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Observation is reported for a structure near the J/y/ij) threshold in B + — > J/yf(j)K + decays pro- 
duced in pp collisions at yfs = 1.96 TeV with a statistical significance of beyond 5 standard 
deviations. There are 19 ±6 events observed for this structure at a mass of 4143.4l^(stat)± 
0.6(syst) MeV/c 2 and a width of 15.3+g°j 4 (stat) ± 2.5(syst) MeV/c 2 , which are consistent with 
the previous measurements reported as evidence of the 7(4140). 
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Figure 1: The mass distribution of J/y/<pK + ; the 
solid line is a fit to the data with a Gaussian signal 
function and linear background function. 
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Figure 2: The mass difference, AM, between 
fl + fl~K + K~ and in the B + mass window is 

shown as the black histogram. The red histogram is 
the AM distribution from the data in the B sideband. 



Recently, evidence has been reported by CDF for a narrow structure near the J / y/0 threshold, 
named 7(4140), in B + — > J /ytyK + decays produced in pp collisions at yfs = 1.96 TeV The 
structure is the first charmonium-like structure decaying into two heavy quarkonium states (cc and 
ss) which is a candidate for exotic mesons In this note, we report an update on a search for 
structures in the J/y<p system produced in exclusive B + — > J/y<pK + decays with J j\\f — > pi + pL~ 
and (p — > K + K . This analysis is based on a sample of pp collision data at yfs = 1.96 TeV with an 
integrated luminosity of about 6.0 fb -1 collected by the CDF II detector at the Tevatron. The CDF 
II detector has been described in detail elsewhere [Q]. In this analysis, J/y— > At + Ai~ events are 
recorded using a dedicated three-level dimuon trigger. 

The invariant mass of J/\j/(f)K + in the current dataset, which includes those used in the pre- 
vious analysis after applying the same requirements used in the previous analysis [jl|], is shown in 
Fig. [I]. A fit with a Gaussian signal function with its rms fixed to the value 5.9 MeV/c 2 obtained 
from Monte Carlo (MC) simulation [Q] and a linear background function to the mass spectrum of 
// Y<j>K + returns a B + signal of 1 15 ± 12(stat) events. For a luminosity increase by a factor of 2.2, 
the yield increase was 1.53, reduced by trigger rate-limitation at higher average luminosity. We 
then select B + signal candidates with a mass within 3a (17.7 MeV/c ) of the nominal B + mass. 
We define those events with a mass within [-9,-6] a or [6,9] a of nominal B mass as B sideband 
events. Fig. || shows the mass difference, AM = m(jj. + iJ.~K + K~) — m(/i + /x _ ), for events in the 
B + mass window as well as in the B sideband in our data sample. The comparison of the AM 
distribution in the B mass window for the dataset used in this analysis (6.0 fb -1 ) and the dataset 
used in the previous analysis (2.7 fb -1 [[I]]) is shown in Figure |3[ 

The same model is used to examine the F(4140) structure as described in reference [|l|]. We 
model the enhancement by an 5-wave relativistic BW function ^ convoluted with a Gaussian 
resolution function with the r.m.s. fixed to 1.7 MeV/c 2 obtained from MC, and use three-body 
phase space [|]] to describe the background shape. Even though we exclude the high mass region to 
avoid the B s contamination, there is still a small contribution in the region of interest. We obtained 
the AM shape from B s contamination and fix the AM shape obtained from B s MC simulation, and 
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Figure 3: The AM distribution in the B mass win- 
dow for the data used in the current analysis (6.0 
fb -1 ) is shown as the black histogram, and the same 
distribution for the data in the previous analysis(2.7 
fb 1 l|j]])i s shown as the red dashed histogram. 
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Figure 4: The mass difference, AM, between 
fl + fl~K + K~ and in the B + mass window is 

shown as a solid black histogram for the data. The 
dotted curve is the predicted three-body phase space 
background contribution, the dash-dotted curve is 
the predicted B s contamination (fixed to 3.3), and 
the solid red curve is the total unbinned fit where 
the signal PDF is an S-wave Breit-Wigner convo- 
luted with the resolution (1.7 MeV/c 2 ). 



the yield to 3.3, scaled from the B s — > J/y<p yield in data. An unbinned likelihood fit to the AM 
distribution, as shown in Fig. ||, returns a yield of 19 ±6 events, a AM of 1046.7+ 2 q MeV/c 2 , and 
a width of 15. 3^ MeV/c 2 . 

We use the same simulation process as in Reference [ffl], based on a pure three-body phase 
space background shape to determine the significance of the 7(4140) structure. We performed a 
total of 84 million simulations and found 19 trials with a a/— 2\n(J£Q / J£ max ) value greater than 
or equal to the value obtained in the data (5.9), as shown in Fig. ||, where J£q and 5£ max are the 
likelihood values for the null hypothesis fit and signal hypothesis fit. The resulting /rvalue is 
2.3 x 10~ 7 , corresponding to a significance of greater than 5.0a. 

The mass of this structure, including systematic uncertainty, is measured as 4143.4+ 2 'o(stat) ± 
0.6(syst) MeV/c 2 after including the world-average J ' j\\f mass. The relative efficiency between B + — >■ 
7(4140)^+, 7(4140) -> J/\j/(j> and B+ -> J/y$K+ is 1.1 assuming 7(4140) as an S-wave struc- 
ture and B + phase space decays. Thus the relative branching fraction between B + — > 7(4140)^+, 
7(4140) -> J/y§ and B + -> J/y$K + including systematics is 0.149 ±0.039(stat) ±0.024(syst). 

An further excess above the three-body phase space background shape appears at approxi- 
mately 1.18 GeV/c 2 in Fig. |l| (b). Since the significance of 7(4140) is greater than 5a, we fit to 
the data assuming two structures at AM of 1.05 and 1.18 GeV/c 2 as shown in Fig. ^. The fit to 
the data with the same requirements as in the previous paper [[J] returns a yield of 20 ± 5 events, a 
AM of 1046.7^ 2 J MeV/c 2 , and a width of 15.0+^ MeV/c 2 for the 7(4140), which are consistent 
with the values returned from a 7(4140)-only signal fit. The fit returns a yield of 22 ± 8 events, a 
AM of 1 177.7+|-* MeV/c 2 , and a width of 32.3+^ ^ MeV/c 2 for the structure around AM of 1.18 
GeV/c 2 . The significance of the second structure, determined by a similar simulation is 3.1a. 
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Figure 5: Distribution of -21n(Jz? / -^max) for 84 
million simulation trials. The p-value obtained from 
counting is 2.3 x 10~ 7 , corresponding to a signifi- 
cance of 5.0(7. 



Figure 6: The mass difference, AM, between 
H + H~K + K~ and n + ji~ , in the B + mass window. 
The dotted curve is the background contribution, the 
dash-dotted curve is the B s contamination, and the 
red solid curve is the total unbinned fit assuming two 
structures. 



In summary, the growing B + — > J/\j/<pK + sample at CDF enables us to observe the 7(4140) 
structure ^ with a significance greater than 5a. Assuming an S- wave relativistic BW, the mass and 
width of this structure, including systematic uncertainties, are measured to be 4143.4+ 2 Q(stat) ± 
0.6(syst) MeV/c 2 and 15.3+g°j 4 (stat) ± 2.5(syst) MeV/c 2 , respectively. The relative branching 
fraction between B + -> Y(4140)K + , 7(4140) -> J/yQ and B + -> J/y(j)K + including systematics 
is 0.149 db 0.039(stat) ± 0.024(syst). We also find evidence at 3.1a level for a second structure 
with a mass of 4274.4+^ (stat) MeV/c 2 , a width of 32.3™ (stat) MeV/c 2 and a yield of 22 ± 8. 
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